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We describe a spin-echo method for ultracold neutrons (UCNs) confined in a precession chamber 
and exposed to a |-Bo| = 1 pT magnetic field. We have demonstrated that the analysis of UCN spin- 
echo resonance signals in combination with knowledge of the ambient magnetic field provides an 
excellent method by which to reconstruct the energy spectrum of a confined ensemble of neutrons. 

The method takes advantage of the relative dephasing of spins arising from a gravitationally induced 
striation of stored UCN of different energies, and also permits an improved determination of the 
vertical magnetic-field gradient with an exceptional accuracy of 1.1 pT/cm. This novel combination 
of a well-known nuclear resonance method and gravitationally induced vertical striation is unique 
in the realm of nuclear and particle physics and should prove to be invaluable for the assessment 
of systematic effects in precision experiments such as searches for an electric dipole moment of the 
neutron or the measurement of the neutron lifetime. 

Keywords: ultracold neutrons, energy spectrum, magnetic field gradient, spin-echo 


Spin-echo resonances, first observed by Hahn in 1950 [1], have proven to be a very powerful tool in nuclear 
magnetic resonance (NMR) experiments for identifying different sources of depolarization in spin-polarized 
samples. In particular, the 7r/2 - tt pulse sequence proposed by Garr and Purcell [2] is widely applied 
to distinguish the intrinsic spin-lattice or spin-spin coherence time T 2 from the global coherence time TJ" 
dominated by dephasing of processing spins due to spatially varying magnetic fields. Another established 
application of the spin-echo technique is in neutron scattering [3] where it is used to resolve correlation times 
below 1 ps. Essentially, the spin-echo technique can be applied in solid and soft matter samples where the 
intrinsic component is long compared to dephasing, which can even be the case for T 2 <C 1 s. Two similar 
components of depolarization [4] exist in experiments using UGNs [5-11] where coherence times of several 
hundred seconds have been achieved. The defining characteristic of UGNs is that they are reflected from 
appropriate material surfaces at all angles of incidence, even at room temperature, conserving kinetic energy 
as they do so. They can therefore be confined for observation times that are commensurate with neutron 
/3-decay: inelastic scattering with phonons from the storage-vessel walls contributes F « 10“®s“* to their 
loss rate. This peculiarity also makes them ideal for neutron lifetime measurements. A measurement of the 
evolution of the energy spectrum in these UGN storage experiments can give insight into energy dependent 
systematic effects and could help to understand the disagreement between beam and storage experiments [12]. 

Ultracold neutrons have kinetic energies Akin of the same order of magnitude as their gravitational potential 
energy Ep^t above the lower confining surface. For any given UGN E = Ekin+Epot is constant during storage, 
although for the ensemble the average energy decreases over time because faster neutrons have a higher loss 
rate. These energy-dependent losses result in a time dependence of both depolarization terms, making T 2 a 
function of the free-precession time. For neutrons mg « 1.03neV/cm, so, e.g., UGN with kinetic energy of 
200 neV can rise about 200 cm in the gravitational field. Therefore, it is convenient (and common) to refer to 
the energy e of UGN in terms of the maximum height (in cm) attainable within Earth’s gravitational field, 
and we do so here during our analysis, although we revert to E (in neV) for our final results. 
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Two channels of field-induced depolarization of UCNs exist: the known intrinsic depolarization due to 
inhomogeneous magnetic fields [13-18], and an energy-dependent relative dephasing for different energies 
in a vertical magnetic-field gradient. The latter shows an interplay between gravitationally defined spatial 
distributions of the neutrons and a vertical magnetic-field gradient dBz/dz and has only recently been 
studied theoretically [4, 19, 24]. It is the subject of this investigation which uses a new form of spin-echo 
spectroscopy that also circumvents effects provoked by the time-dependent softening of the UCN spectrum. 
Prior to this, the spin-echo resonance method was used in a study of Berry’s phase with UCNs to cancel any 
dynamic phase shift [20] with a coherence time of approximately 850 ms. 

We have measured UCN spin-echo signals using an apparatus dedicated to the neutron electric dipole 
moment (nEDM) experimental] with a vertical magnetic field of Bq « (0,0,—l)pT and depolarization 
times of more than 450 s. These measurements, together with first-order knowledge of the magnetic-field 
gradients, were used to determine the energy spectrum of the UCN after storage as well as a precise value 
for a common offset to the vertical magnetic-field gradients. The technique has the potential to analyze and 
correct for the so-called geometric-phase effect [15, 35], which is currently the most important systematic 
effect in experiments searching for an nEDM. 

The equilibrium density distribution of UCNs as a function of height h above the lower confinement surface 
in a bottle with diffusely scattering walls is p{h) = / 9 ( 0 ) • \/l — h/e [22]. This leads to a center-of-mass offset 
of H/2 — {h{e)) relative to the center plane of a cylindrical precession chamber of radius r and height H, 
where 
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is the time-averaged height of a neutron [4] with 77 = (1 — 

In the presence of a vertical magnetic-field gradient, neutrons with an energy e will precess with an average 
Larmor frequency of 


^(e) = 7n (± {Bo) - ^ (i?/2-(Me)))) , (2) 

where ± {Bq) is the volume-averaged magnetic field magnitude in the precession chamber for Bq up (-I-) or 
down (—), and 7 n = —27r-29.1646943(69) MHz/T [25] is the gyromagnetic ratio of the neutron. This induces 
an energy-dependent relative phase <I)(e) = w(e)T, where T is the free spin-precession (ESP) time. It leads to 
a relative dephasing which, when averaged over the energy spectrum, can be interpreted as a gravitationally 
enhanced depolarization. This relative dephasing can be studied by applying a7r/2-7r-7r/2 spin-echo 
sequence, varying the time ti at which the 7r-pulse of duration is applied while keeping constant the total 
duration T + between the end of the first and start of the second 7r/2 pulse. This second 7r/2 pulse, at the 
end of the sequence, is required to transfer information about the spin-precession phase onto the measurable 
longitudinal spin component. For convenience, we analyze the system in a rotating frame precessing with 
Wo = In'{Bo) and Bo>0 inside the storage volume, and we only consider a linear gradient dB^fdz along 
the primary Bq field axis z. Therefore, in the rotating frame, the spins of UCNs with energy E precess with 
average frequency (to first order in dB^jdz) 


Wr(e) = -7„^-(i?/2-(Me)))- (3) 

Figure 1 illustrates the evolution of the spins in the rotating frame. The resonance condition under which 
the spins of all different energies refocus is established when T = 2ti. At this time we measure the highest 
polarization, which is essentially the polarization {a{T, e))e without the gravitationally induced dephasing. 
If we change the start time ti of the 7r-pulse to values shorter than r/2, the resonance condition will have 
been met before we apply the second 7r/2-pulse. As the spins of the different energy classes then proceed to 
fan out again, the final polarization decreases. Similarly, if the 7r-pulse is applied at times larger than T/2 
the refocusing is not yet complete by the time the second 7r/2-pulse is applied. 

For a particular UCN energy e and fixed T the observed polarization after the second 7r/2-pulse can be 
written as 


P{e,ti) = a{T,e) ■ cos (wr(e) (T - 2ti)). 


( 4 ) 
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Figure 1: Illustration of a UCN spin-echo measurement of duration T = 2ti+t{ in the frame rotating at frequency 
cJrf = u!o for a negative gradient dBzjdz: (a) a polarized ensemble of UCNs is loaded into the cell, (a-b) an initial 
7r/2-pulse tips all spins (black arrow) into the equatorial plane, (6-c) the UCN spins precess with uJy{E) for ti, and fan 
out. Low-energy UCNs (blue) see a larger field for a negative gradient dBzjdz than higher-energy UCNs (magenta), 
whereas spins that precess at tuo (i.e. with no center-of-mass offset) are stationary and are oriented along y. A 7r-pulse 
(c-d) then flips the spins around the a;-axis, after which (d-e) the UCNs continue to precess in the same direction, 
eventually refocusing at 2t\ (black arrow). As they continue to precess beyond this time they fan out again (e-/) 
until finally {f-g) a second 7r/2 pulse is applied. All pulses used for spin manipulation have the same field strength 
Bx applied along the same axis {x). A single oscillator is used, with the output being gated on only during the blue 
time windows; hence, all of the pulses are phase coherent with each other. All energy-dependent effects during the 
pulses cancel, since a full rotation about 27r is made and only dephasing during FSP has to be taken into account. 


By integrating over the UCN energy spectrum p(e, T) at the time of polarization analysis T the observed 
polarization is then 


P{ti) = y'a(T,e)-cos [a;r(e) (T-2ti)]p(e,T)de, (5) 

where p(e, T) is normalized to 1 for a fixed T. 

We have used the MCUCN package [23] with a detailed model of the nEDM experiment at Paul Scherrer 
Institute (PSI) [21] {H — 12 cm, r = 23.5 cm) and simulated the measurement with a 7r/2 - tt - 7r/2 spin-flip 
sequence (tip - flip - tip) for a fully polarized initial population of UCNs. This simulation was carried out 
for two different UCN energy spectra p(e), parametrized as in equation (7) below, each with three different 
constant vertical gradients of dB^/dz = 100, 200, and 400pT/cm for a free precession time of T = 92 s. 
The resulting simulated spin-echo resonances are shown in the inserts of Fig. 2 which can be described by 
equation (5). In particular we used a{T, e) = ap sxp (—T • r 2 (e)), with a transverse relaxation rate as derived 
in [24] of 
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where v{e), and ’H(e) = e V e < and (e) = H \/ e > H are the energy-dependent mean velocity, and the 
effective height of the UCN as introduced in [22] for a cylindrical cell; a is a free fit parameter. We found it 
useful to model the UCN spectrum for a fixed precession time T as 
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Figure 2: Simulated spin-echo signals and extracted spectra. The inserts show the simulated spin-echo signals using 
MCUCN[23]. On each plot three resonances are shown, using the gray input spectra with three different vertical field 
gradients: blue lOOpT/cm; green 200pT/cm; and red 400pT/cm. The points with their uncertainties come from the 
simulation. The red lines on the main plots are the best simultaneous fits to the simulated spectra. The colors (from 
bright yellow to dark green) represent the value for each line; dark green corresponds to a 68.3% C.L while bright 
yellow lines have x^ values close to the minimum. The original input spectra are represented by the gray areas, which 
give center-of-mass offsets of 1.93mm (Case A) and 6.46mm (Case B). 


A- e^/2 / I I 

l + e(^o-e)/b^ ' [l + e(e-(Eo+Ei))/c^ I ^ g(E2-Eo-Ei)/c=‘ 

where i? 0 ) b and Eq + Ei, c are the points of inflection and slopes of the rising and falling energy edges 
respectively, while Eq + Ei + E 2 is the high-energy cut-off. A is chosen such that /g^^p(r, e)de = 1. 
This energy spectrum is based on a very general distribution n{e)de oc e^^^de from the low-energy tail of 
a Maxwell-Boltzmann distribution, allowing for low- and high-energy cut-offs. A simultaneous nonlinear 
least-squares fit (algorithm: ‘trust-region-reflective’) to the signals of all three gradients yielded a, 6, c, Eq, 
El, and E 2 for each of the two spectra. The upper and lower confidence intervals were estimated by carrying 
out random scans of the parameter space around the best value with < Xumit 5 x 10® samples; yellow 
to dark green curves in Fig. 2): 



xLit = xLin (^1 + n-p)^ , (8) 

where Fa is the inverse cumulative function of the statistical F-distribution for an a = 68.3% confidence 
level with p = 6 free parameters and (n — p) = 143 degrees of freedom [26]. For the numerical integration 
of equation (5) we used energy bins of equal height: Ah{e) = 0.02 mm, which is sufficiently small that any 
systematic error from the finite energy bins would be less than the statistical resolution of each data point. 
Figure 2 compares spectra extracted from the fit to the input spectra of the simulation. We note that the 
results become less reliable at higher energies; this is to be expected, since once the neutrons populate the 
bottle more or less uniformly it becomes increasingly difficult to discriminate between them. 

We measured three spin-echo profiles (see Fig. 3) using the apparatus described in Ref. [27], improved 
with an array of 16 cesium-vapor magnetometers (CsM) [28] and a simultaneous spin-analyzing system [29]. 
The magnetic-field measurement with the CsM array permitted us to calculate the magnetic-field gradients 
using a polynomial held decomposition up to second order, as described in Ref. [31]. Data were taken for 
a constant FSP duration of T = 216s, with three different trim-coil current settings for Bq < O(down), 
thereby applying vertical gradients of dB^/dz = —18, 10, and 38pT/cm. The uncertainties of the gradient 
measurements are dominated by systematics due to individual sensor offsets, and amount to 12-14 pT/cm 
(see [31, 32] for a determination of these gradient errors). However, under our conditions, such individual but 
constant offsets lead to a constant offset Gz of the true dBz/dz that can be determined to within l-2pT/cm 
using the treatment described below. 
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Figure 3: Measured UCN spin-echo signal for three nominal magnetic field gradients. Red squares: dB^/dz = 
— 18pT/cm; black diamonds: lOpT/cm; blue circles: 38pT/cm. The lines represent the best fit to the data of 
equation (5), using (7) and (9). The residuals are r = P{ti)m — m) exp ■ 


The measured initial polarization ao = 0.861(2), 0.865(2), 0.861(2) was prepared by passing UCNs from 
the UCN source at the PSI [30] through a 5T solenoid magnet. After a filling time of 24s, the UCN shutter 
in the ground electrode was closed. For each setting, one fixed radio frequency Wrf = 30.2075,30.2031, 
and 30.2013 Hz was used. The ^®®Hg co-magnetometer measured the mean magnetic field Bq via WHg = 
7Hg-(i?o) [33]. 

A full phase coherent tt/2 - tt - tt/2 spin-flip sequence was applied to the neutrons in the cell before the 
UCN shutter was then opened again, and the UCN were detected in a pair of spin-state-specific detector 
channels. A measurement consisted of several cycles with ti increasing in steps of 9 s. Figure 3 shows the 
three measured spin-echo resonances. As for the simulated cases, equation (5) was fitted to the data, using 
the same fitting procedure and in addition randomly varying the initial start values. The frequency lo^ from 
equation (3) was corrected for the difference Auj = uj^f — |7n|/7Hg<^Hg and was also extended by a term for 
the second-order gradient in z, which was available from the field measurement: 


Wr(e) =7n 


■ {H/2 - (/i(e))) -I- 


dz^ 


.{{h^)-{{H/2-h{e)f)) 


— Auj, 


(9) 


where (h^) = z^/Hdz = 12 cm^ is the spectrum-independent quadratic expectation value. The fit 
parameter Gz was included to accommodate the above-mentioned common gradient offset. As the ob¬ 
served polarization at = T/2 was the same for all three measurements, we used an averaged com¬ 
mon F 2 (e) for the fit. Figured shows the extracted energy spectrum. The best fit parameters were: 
Gz = 5.8(1.l)pT/cm, a = 10.0(4), b = lt°;VeVi/2, c = 2.7^1 neV^/^, Eq = 7.7VjneV, and Ei = 
28.7l'l4'with E 2 fixed at 200 neV. Lower bounds for b = IneV^/^ and c = 2.5 neV^/^ were chosen 
to give smooth physical edges while remained essentially constant for even lower values. The average 
energy (E) of the measured spectrum is 27.6^V neV, corresponding to a mean velocity at the base of the con¬ 
tainment vessel of T = 230^:)® cm/s. The observed intrinsic depolarization time at U = T/2 is T 2 = 506(26) s, 
while the calculated dephasing times, correcting for the frequency offset Aw are = 449(8), 439(10), and 
279(9) s for the measurements at nominal dBz/dz = —18, 10 , 38pT/cm respectively. Note that a = 10.0(4) 
indicates that the intrinsic depolarization is a factor 10 higher than calculated which we attribute to higher 
order magnetic-field gradients or magnetic impurities of the walls of the precession cell. 

The extracted spectrum-averaged UCN to i®®Hg center-of-mass offset hos = ® mm is larger using this 

method than the value 2.35(5) mm [31] deduced from a measurement of the ratio of the precession frequencies 
of i®®Hg atoms to neutrons, R = Wn/wHg, after 180s of storage (see Ref. [31]). This difference is explained 
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Figure 4: Normalized energy spectra of stored UCNs after 220 s storage. The red line displays the energy spectrum 
obtained from the best global fit to equation (5), using (7) and (9), to all three data sets, compare Fig. 3. The 
additional spectra shown are those having the highest and lowest mean energy (brown) and the highest and lowest 
center-of-mass offset (blue) within a 68.3% confidence interval. 


by the effect on the measurement of R of gravitationally enhanced depolarization, as discussed in [4, 24], as 
well as a softening of the spectrum for longer storage times. 


We have demonstrated, using simulations, that the analysis of UCN spin-echo resonance signals in combi¬ 
nation with knowledge of the ambient magnetic field provides an excellent method by which to reconstruct 
the energy spectrum of a confined ensemble of neutrons. The method takes advantage of the relative de¬ 
phasing of spins arising from a gravitationally induced striation of stored UCN of different energies. Based 
upon this, measured UCN spin-echo signals have been analyzed. They were taken in a controlled magnetic 
environment which allowed for long inherent coherence times and permitted a FSP time of 216 s. From these 
measurements it was possible not only to extract the energy spectrum of stored UCNs, but also to determine 
a common gradient offset with a resolution of l.lpT/cm, and thus to determine the magnetic-field gradient 
at this level of accuracy. We plan to improve the energy resolution at higher energies (50 to 200 neV) by 
repeating this measurement with a larger bottle of ~ 40 cm. Using an array of vector magnetometers, 
as described in e.g. Ref. [34], will further improve this technique. Obtaining the optimal resolution will also 
require an improved understanding of the intrinsic-depolarization processes a(T,e), particularly in the limit 
of very low energies, i.e. E <H, including effects of different specularities of wall reflections and higher order 
magnetic-field gradients. 

This technique will improve the estimation of energy-spectrum dependent systematic effects in high- 
precision experiments such as searches for the neutron electric dipole moment [21] and for spin-dependent 
forces [32]. Measurements for different T will give access to the evolution of the spectrum p{e,T) during 
storage, which may help to solve the neutron lifetime controversy [12]. Alternating measurements with 
and without the 7r-pulse at U = T/2 could also provide a powerful method to improve the magnetic-field 
homogeneity of typical UCN spin-precession experiments, such as, for example, neutron EDM searches, 
where it is known that magnetic field gradients can cause false EDM signals when a co-magnetometer is 
used [15, 35]. Using the gradient-offset extraction from a spin-echo measurement in combination with a 
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time-resolved gradient measurement by the CsM array may allow one to correct directly for these false-EDM 
signals in the future. 
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